Pairing and radio-frequency spectroscopy in two-dimensional Fermi gases 
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We theoretically study the normal phase properties of strongly interacting two-component Fermi 
gases in two spatial dimensions. In the limit of weak attraction, we find that the gas can be 
described in terms of effective polarons. As the attraction between fermions increases, we find a 
crossover from a gas of non-interacting polarons to a pseudogap regime. We investigate how this 
crossover is manifested in the radio-frequency (rf) spectroscopy. Our findings qualitatively explain 
the differences in the recent rf spectroscopy measurements of two-dimensional Fermi gases [Sommer 
et al, Phys. Rev. Lett. 108, 045302 (2012) and Zhang et al., Phys. Rev. Lett. 108, 235302 (2012)]. 
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I. INTRODUCTION 

Pairing of fermions in two spatial dimensions has 
become one of the main themes in condensed-matter 
physics due to the discovery of high-temperature su- 
perconductivity in copper oxide compounds (l|. Ultra- 
cold gases of fermionic atoms provide a flexible platform 
for testing the various pairing mechanisms. Advances 
in manipulating atoms in optical lattices can ultimately 
lead to direct simulation of materials relevant to high- 
temperature superconductivity. To achieve this mile- 
stone, one first needs to understand pairing in simpler 
systems, such as Fermi gases confined to continuous two- 
dimensional (2D) geometries. Several recent experiments 
have probed different aspects of strongly interacting 2D 
Fermi gases, ranging from polaron physics 2hJ] to pair- 
ing and confinement-induced molecules [5H8J- Another 
interesting direction is a dimensional crossover from two 
dimensions to three dimensions [5|, |6[ • 

One of the intriguing aspects of both high-temperature 
superconductors and atomic Fermi gases is the possible 
existence of preformed Cooper pairs above the super- 
fluid transition temperature [l|, I9l4l4|. This regime is 
referred to as the pseudogap phase, and in the context of 
Fermi gases, the properties of the pseudogap regime have 
been experimentally probed in both two-dimensional and 
three-dimensional (3D) systems [71 [R* - [l7| . The nature of 
the normal state of a strongly interacting Fermi gas is still 
an open question since it can also be interpreted in terms 
of the Fermi-liquid framework [l8j. Furthermore, recent 
experiments have suggested that a 2D Fermi gas can be 
described as a gas of noninteracting polarons even in the 
absence of any population imbalance [J] , while other ex- 
periments indicate the existence of confinement-induced 
molecules in the normal phase [f|. 

In this work, we probe the properties of the normal 
state of strongly interacting 2D Fermi gases. In partic- 
ular, we aim to qualitatively explain the differences in 
the recent experiments performed in the Zwierleingroup 
at Massachusetts Institute of Technology (MIT) Q and 
in the Thomas group at North Carolina State University 
(NCSU) Q . We find that the normal state of weakly 



attractive Fermi gases can be described in terms of effec- 
tive polarons, whereas strongly attractive Fermi gases are 
characterized by a Bardeen-Cooper-Schrieffer (BCS) -like 
effective dispersion and suppression of a single-particle 
density of states (DOS) near the Fermi energy. Such 
characteristics are commonly attributed to the pseudo- 
gap regime Q3 E E GJ, EE H Hfll . Our analysis sug- 
gests a crossover from a gas of non-interacting polarons 
to the pseudogap regime with increasing attraction and 
qualitatively explains the different experimental results 
reported in Refs. @, @ • We note that the earlier experi- 
ments reported in Ref. [2( regarding molecule formation 
in 2D Fermi gases have already been interpreted in terms 
of dynamically created polarons [21( and therefore we do 
not discuss this experiment in the present work. 

The experiments reported in Refs. [J, Q used radio- 
frequency (rf) spectroscopy to probe the properties of 
2D Fermi gases. Both experiments can be schemati- 
cally described by an interacting, population balanced 
initial mixture of spin-f and spin- 1 atoms. A short rf 
pulse is applied to convert spin- J, atoms into a final-state 
|/) and the number of converted atoms (or, equivalcntly, 
atoms remaining in state \ D) is subsequently measured. 
By varying the detuning of the rf pulse from the bare 
atomic transition, one obtains information on the single- 
particle excitation spectrum. For this reason, rf spec- 
troscopy has been extensively used to probe the pai ring 
mechanisms in strongly interacting Fermi gases [171 122| - 
|26| . Although both experiments probed, in principle, the 
same initial system, the final conclusions were quite dif- 
ferent. The MIT experiment [6[ was interpreted in terms 
of confinement-enhanced pairing [27|, whereas the NCSU 
experiment J4| found the resonances in the rf spectra to 
be best explained by transitions between polarons in the 
initial and final-states. 

To analyze the different contributions to the experi- 
mentally measured rf spectra, one needs to take into ac- 
count interactions between the initial-state atoms (| f) 
and | \)) and the final-state atoms (| f) and |/)). On 
the other hand, interactions between atoms in states | \) 
and \f) are irrelevant since the rf pulse coherently ro- 
tates the atomic spin [28|. We note that both the MIT 



and the NCSU experiments correspond to relatively weak 
final-state interactions [J, Q ■ Therefore we consider the 
final-state interactions only phcnomenologically and con- 
centrate on the strong initial-state interactions. The con- 
tribution of the final-state interactions to the MIT exper- 
iment has recently been discussed in Ref. [20|. In Sec.lPH 
we describe the non-self-consistent T-matrix formalism 
which we use to calculate the spectral properties of the 
initial- and the final-state atoms. Section Mil describes 
the different schemes for calculating the chemical poten- 
tial of the population balanced initial state. The quasi- 
particle spectrum of the initial- and the final-state atoms 
is discussed in Sec. HV] The rf spectroscopy correspond- 
ing to the experiments in Refs. 0,0] is studied in SeclVl 



and concluding remarks are presented in Sec. [VT 



II. T-MATRIX AND THE LADDER 
APPROXIMATION 

A. Initial state 

We study the properties of the population balanced 
initial state using a non-self-consistcnt T-matrix approx- 
imation where the T-matrix is computed by summing 
over all ladder diagrams. This approximation has been 
extensively utilized in the literature to study the proper- 
ties of Fermi gases in the normal phase jlCX IllL ll3L [30M33| . 
Alternative theoretical approaches have been discussed in 
Ref. [14| ■ The dressed Green's function can be calculated 
from the Dyson's equation, 

Q~ 1 {k,iu n ) = G^(k,iuj n ) - T,a.(k,w n ), (1) 



where the bare 
Q^^k.iuJn) = icu n - 



Green's function is given by 
s k + /i CT , and a =t, 4-- The fermionic 
Matsubara frequencies are defined as ui n = (2n + l)ir//3 
and j3 = l/k B T. Furthermore, we have set H = 1. 
Within the non-self-consistent ladder approxima- 
tion [10|, [U Oil [301 - 1331 ] , the self-energy is given by the 
T-matrix and the bare Green's function, 



£ CT (fc,iwJ 



-T 



dq 



(2tt) 2 
0o.- 



%(<L^m)x 



fe, iCl, 



"O- ( 2 ) 



Notation — a indicates a spin opposite of a and fl n = 
27m//3 denote the bosonic Matsubara frequencies. The 
many-body T-matrix can be calculated from the Bethe- 
Salpeter equation using the ladder approximation and 
the bare Green's function 

t 77. o S = 7T + Xo(g,*"m), (3) 

such that the polarization operator is given by 
1 ^ f dk 



x (q,ittr 
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(2tt)2 
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x Go,i(k,iu) n ). (4) 



The ladder approximation for computing the self-energy 
is illustrated in Fig. [TJ 

We assume the gas has equal densities of spin- 4- and 
spin-f particles, which implies /ij, = //-(- = pb. To can- 
cel the ultraviolet (UV) divergence associated with the 
polarization operator, we use the vacuum T-matrix 



r (w) = 



47r/m 



ln(e B /w) + in 



(5) 



where e B = 1/ma^ > is the two-body bound-state 
energy. The many-body T-matrix takes the form 
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where the regularized polarization operator XQ° s (q,iJ7 m ) 
is given by 
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and n F {ui) denotes the Fermi function. In this work, we 
are interested in finite temperatures above the supcrfluid 
phase-transition temperature, and the polarization oper- 
ator in Eq. ([7]) has to be computed numerically. 
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FIG. 1. Ladder approximation for the self-energy. We restrict 
the analysis to the fluctuation part since the exchange part 
is important only in the case of magnetic ordering. Further- 
more, we assume equal densities of spin- J, and spin-'f atoms, 
which makes the diagrams interchangeable with respect to 
spin-t and spin- J,. 

After the analytic continuation iio n — > uj + z0 + , the re- 
tarded self-energy E, (qr, uj + i0 + ) can be computed from 
Eq. ^ using contour integration. Since the T-matrix 
may have poles away from the real axis, we cannot di- 
rectly apply the spectral representation of the T-matrix 
to compute the self-energy. However, assuming that the 
T-matrix has poles either on the real axis or symmetri- 
cally with respect to the real axis, we arrive at a simple 



result for the imaginary part of the self-energy, 



polarization operator is given by 



/dfc 
— — [n B (u + £ k - m) + n F {e k - //)] 

x1mT n (q + k,u + £ k -iJ,), (8) 

where n B (uj) denotes the Bose distribution. We have 
numerically verified that, for the parameter regime con- 
sidered in this work, the T-matrix has poles at most on 
the real axis. The real part of the self-energy is calculated 
using the Kramers-Kronig relation 



r°° dz 

ReZ i (q,u J )=V / - 



dz ImS,(qr,z) 



(9) 



where V denotes principal- value integration. 

In 3D systems the appearance of poles in the T-matrix 
with a nonzero imaginary part is associated with the on- 
set of the superfluid phase. This is the Thouless crite- 
rion [3J], and it can be used to identify the superfluid 
transition temperature T c . In 2D Fermi gases, the super- 
fluid transition takes place via the Bcrczinskii-Kostcrlitz- 
Thoulcss (BKT) transition, which is not captured by our 
non-self-consistent T-matrix theory [9| . In order to access 
the BKT physics, one needs to either consider the vertex 
corrections to the T-matrix formalism |9| or explicitly in- 
clude phase fluctuations to a mean- field formalism |35J |. 
Our non-self-consistent T-matrix theory should therefore 
be used at temperatures above T c where it provides a rea- 
sonable description of the system. 

The retarded Green's function is obtained by analytic 



continuation Gf-(q, 



UJ) 



QAq,iu r , 



i0 + ) from 



the thermal Green's function given by the Dyson's equa- 
tion ([TJ). The corresponding spectral function is defined 
as AJq, u) = —2 Im Gf~(q, ui), and it satisfies a sum rule 

I^oo fir Al(9> w ) = 1- We use this sum rule to verify 
the consistency of the numerical calculation. To obtain a 
physically correct spectral function for a strongly inter- 
acting initial gas of spin-'f and spin- J, atoms, the chemi- 
cal potential /x has to be determined self-consistently. In 
the next section, we discuss the different self-consistent 
schemes for computing \x. 



B. Final state 

In Rcf . y , the experimental findings were interpreted 
in terms of transitions between initial- and final-state po- 
larons. In order to test this hypothesis, we consider the 
properties of the final-state atoms when they are dressed 
by the interactions with the initial state atoms. We take 
the interactions between atoms in hyperfine spin states 
| f) and | /) into account, using again the T-matrix for- 
malism [21j. For the final-state atoms, the regularized 



Xo eg (<?,iO m ) = 

dfc n F (e k+q/2 -^)+n F (e k _ q/2 + A ff ) 
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zt fc 2 t q 
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21|, and 

with 



where Aff is the level splitting between states | 'f) and 
|/). We have also denoted ix = fj,±. Since the final-state 
is initially empty, we have taken the final-state chemical 
potential to be fi f = 0. We use the vacuum T-matrix to 
regularize the UV divergence. This introduces an addi- 
tional interaction parameter e' B characterizing the final- 
state interactions. 

Calculation of the dressed Green's function is analo- 
gous to the polaron problem considered in Rcf. 
we obtain the self-energy E* by convoluting t/o,t 
the many-body T-matrix T^f. Since the T-matrix cor- 
responding to the final-state interactions does not have 
poles away from the real axis 2l| , we can directly a pply 
spectral representation for the T-matrix and obtain [36| 

/dfc 
——2 [n B (u + e k - n) + n F (e k - /j,)] 

x ImT tf (q + k,uj + e k -fi). (II) 

Using Eq. (fTTj) . we calculate the real part using the 
Kramers-Kronig relation ©. In Rcf. |21(, the self-energy 
of the final-state atoms was calculated at T = without 
invoking the Kramers-Kronig relation. As a check for the 
numerical calculation, we have verified that we reproduce 
the results of Ref. [21j in the appropriate limit. 



III. CALCULATION OF THE CHEMICAL 
POTENTIAL 



In this section, we compare two different schemes 
for calculating the chemical potential of the initial-state 
atoms. The first one is the Nozieres-Schmitt-Rink (NSR) 
approximation J30l - l32l . l37| and the second method is 
based on the number density given by the dressed Green's 
function [1J, [38§ . The motivation for studying the two 
approximations is the fact that the NSR scheme is numer- 
ically much more affordable but expected to be accurate 
only when self-energy corrections are small [38|]. We are 
not aware of any explicit comparison of the two schemes 
for 2D systems. 

Both aforesaid methods start from the thermodynam- 
ical potential f2[E, Q] [3g, where E and Q are self-energy 
and dressed Green's function, respectively. The density 
is given by n = — ^ , and the NSR approximation is ob- 
tained by replacing f2[E,£] — > 51[E,0j_. A more rigorous 
approximation suggested by Serene [38| corresponds to 
il[E, Q] — > f2[E, Gq\- The NSR approximation results in a 



number equation of the form [3ll . |32| . |38|, |3£ 



1 J^ n ^ £ i'^ + 
d f dq 



djj,J (2tt) 5 



+ 00 A. . 

™n B {u)5{q,u), (12) 



where 5(q,uj) is the phase shift of the T-matrix, that is, 
T n (q,u)) = \T n (q,w)\e iS ( q ' u) . The total density can be 
written in terms of the Fermi energy as n = ms F /ir. 
Numerical solution of Eq. (jT2"j) determines the chemical 
potential /i in terms of the Fermi energy e F ~ k F /2m. 
We denote the temperature corresponding to Fermi en- 
ergy by T F . 

A more rigorous alternative to the NSR number equa- 
tion is a full number equation where the density is given 
by the dressed Green's function [ljj, [Tj, [38( 



-Y 

fit? 



-^_g(q, lU ; n )e^\ (13) 



such that Q(q,iui n ) is obtained from the Dyson's equa- 
tion (JT|). Using the spectral representation, we can cast 
Eq. ([HI) into 



n = 2 



(2tt)' 



dw 

2^ 



A(q,u)n F (u). 



(14) 



The NSR number equation is obtained from Eq. ([TB")) if 
Q(q,iu} n ) is replaced with |13l.[38j 

g NSK {q,iuJ n ) = g a (q,iuj n )+g (q,iuJ n )Z(q,iu; n )g {q,iuJ n ). 

(15) 
The NSR approximation is therefore formally reliable 
when corrections due to nonzero self-energy are small. 



fa. 




Before analyzing the 2D case, we note that in 3D (at 
the critical temperature), the NSR approximation and 
the full number equation yield chemical potentials that 
are almost identical JT3J . In Fig. [2] we compare the results 
of the NSR approximation and the full number equation 
for a 2D Fermi gas. We present the results in terms of a 
dimensionless interaction parameter [7|, |2l| , 

V = -^ln(e B /2e F ) = ln(k F a 2D ). (16) 

The 2D unitarity is defined as rj = since perturba- 
tive expansions in terms of I/77 diverge at this point [40| . 
A BCS-type supcrfluid and a Bose-Einstein conden- 
sate (BEC) of tightly bound molecules correspond to 
?/ ^> 1 and ?/ <C —1, respectively [23, I4JJ. We observe 
that the difference between the NSR approximation and 
the full number equation is surprisingly small and the 
NSR theory is in fact a good approximation for the full 
number equation in the regime of parameters relevant to 
this work. As one would expect, the NSR approximation 
works best in the weakly interacting regime where ?/ > 1 
(BCS limit). We will use, however, the chemical poten- 
tial computed from the full number equation (fT4")l in the 
subsequent calculations. In Fig. [31 we show the chemical 
potential as a function of the interaction parameter r\ at 
different temperatures. 




FIG. 3. (Color online) Chemical potential at different tem- 
peratures as a function of the interaction parameter rj. From 
top to bottom: T/T F = 1.0, 1.5, 2.0, 2.5, and 3.0. 



IV. SPECTRAL FUNCTIONS AND 
QUASIPARTICLE SPECTRUM 



FIG. 2. (Color online) Chemical potential as a function of 
temperature for different values of the interaction parameter 
77. From top to bottom: 77 = 1.5, 0.5, 0.25, 0, and -0.25. The 
solid line corresponds to the full number equation [Eq. (|14|l ] 
and the dashed line is the NSR approximation [Eq. 1)120 ]. 



We consider two separate cases of strongly interacting 
Fermi atoms: (A) a balanced mixture of spin-J, and spin- 
t atoms and (B) an impurity problem where an atom in 
a state \f) is immersed in a Fermi sea of spin-f atoms. In 
the context of rf spectroscopy [j,[a|, case (A) corresponds 
to the initial state of the system and case (B) describes 



the final-state atoms interacting with spin-f atoms. In 
case (B) we assume that the impurity only interacts with 
spin-f atoms since the rf pulse transferring atoms from 
state | 4-) to the final-state |/) simply rotates the atomic 
spin. 



A. Initial balanced mixture 

The spectral function A,(k,u>) for different values of 
the interaction parameter rj is shown in Fig.|3]for T/T F = 
0.5 and 1.0. In order to reach qualitative agreement with 
experimentally measured spectral functions [7|, one has 
to average the spectral function over the inhomogeneous 
density to account for the free atoms residing at the edge 
of the trap. Such atoms typically give rise to a peak at 
u = -/i [42l - l44l ] (see also Sec. IVCjl . To understand the 
properties of the resulting rf spectra, we analyze next 
the different contributions to the quasiparticle spectra. 
For 77 = —0.5 and 0, we find two distinct features in 
the spectral function: a broad incoherent band and a 
narrow coherent band. This two-band structure is a first 
indication of the pseudogap regime [RJ, H2H151 \YA Il9l . [20j 
since it suggests that removing a spin-], particle from the 
system corresponds to creating and annihilating a pair of 
quasiparticles as in the usual BCS theory [13J . 

We observe that the threshold for the incoherent 
branch in the spectral function (for 77 = —0.5 and 0) 
is approximately given by 



y th 



(9) 



r(?)+M) 



(17) 



where the bound-state energy £ pa i r (<7) is given by the pole 
of the many-body T-matrix. We find that the poles have 
only a vanishingly small imaginary part. The threshold 
o; th (q) is obtained from a simple argument: to create a 
spin- 4- quasiparticle excitation with momentum g, one 
can create a fermion pair with momentum q + k and re- 
move a spin-f atom with momentum fc. This requires 
energy E = e pair (fc + q) — (e k — m) and the threshold 
is obtained for |fc| ~ 0. In a balanced system, this re- 
sults in an energy threshold w th (q) = £ pa i r (<?) + H for the 
incoherent part of the quasiparticle spectrum. We note 
that a similar argument has been put forward earlier in 
Ref. [H. 

To gain more insight into the two branches of the spec- 
tral function for ri = —0.5 and 0, we fit a BCS-like dis- 
persion [TSEMllil, 



'BCS 



(fc) 



(s k + U-^ + Al g 



(18) 



to the numerical data. The pseudogap in Eq. (fT5]l is de- 
noted by A and U corresponds to a Hartrec shift. The 
spectral function in Fig. 2] is shown relative to the chemi- 
cal potential /j, and therefore the two bands are described 
by an effective dispersion, E±(k) = fi±e BCS (k) (19l.|44J|. 
The spectral function itself is peaked around the effec- 
tive BCS-like dispersion given by e±(k) — ±e BCS (fc). 
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FIG. 4. (Color online) Spectral function in the energy- 
momentum plane. Left panels: T/T F = 0.5. Right panels: 
T/T F — 1.0. The color scale is given by log 10 AAk,uj). For 
rj = —0.5 and r\ — 0, the dashed lines indicate a BCS-like dis- 
persion given by Eq. (JTHJ and the solid line is the threshold 
given by Eq. (|17|l . For 77 = 1.0, the dashed line corresponds 
to the free particle dispersion e h — k 2 /2m and the solid line 
indicates the dispersion of the attractive polaron. 



Results from a least squares fit are shown in Table U and 
we observe that the pseudogap A decreases with in- 
creasing temperature and 77. In general, we find that the 
coherent particle branch tends to be more accurately de- 
scribed by the BCS-like dispersion than the hole branch. 
From Fig. 2] and Table U we conclude that for 77 = —0.5 
and 77 = 0, the pseudogap regime extends at least up to 
T = T F if we define it to correspond to a regime where 
A pg is nonzero. We note that the true pairing gap corre- 
sponding to the supcrfluid phase of a 2D Fermi gas has 
been calculated at T = in Ref. [4l| . 

In the context of 3D Fermi gases, the pseudogap 
regime has been investigated by analyzing the density of 
states [13, [201 ■ For 2D systems, the corresponding DOS 
is given by 



p(u) 



dfc 



A(k,i 



(19) 



In Fig. [SI we show the DOS corresponding to Fig. 2) The 
density of states is measured with respect to the corre- 
sponding density of states of a nonintcracting Fermi gas. 
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FIG. 5. (Color online) Density of states p(u) for (a) r\ = —0.5, (b) 77 = 0, and (c) 77 = 1.0. The upper (lower) curves with 
respect to the left peak correspond to T/T F = 0.5 (T/T F = 1.0). The density of states is measured in the units of density of 
states for an ideal Fermi gas, p = m/h . The parameters are the same as in Fig. [4] 
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TABLE I. Pseudogap A pg and Hartree shift U arising from 
a least-squares fit of a BCS-like dispersion to the numerical 
data in Fig. Q] 



In two spatial dimensions, ideal-gas DOS is a constant, 
p = m/h 2 . In terms of DOS, the pseudogap regime can 
be identified as a regime where DOS becomes strongly 
suppressed at zero energy (with respect to chemical po- 
tential) [13, |20]. Using this criterion, we observe that 
the system is in the pseudogap regime for 77 = —0.5 and 
77 = at least up to temperature T = T F ; sec Fig. [5) On 
the other hand, the pseudogap practically vanishes for 
77 = 1.0 at T/T F = 0.5 and 1.0. Another characteristic of 
the pseudogap regime is "backbending" of the quasiparti- 
cle dispersion in the single- par ticle spectral function near 
the Fermi wave vector k F [lj|. From Fig. HJ we observe 
that the backbending near k F is clearly manifested in the 
spectral function for 77 = —0.5 and 77 = at T/T F = 0.5. 
We note that the backbending in the normal state of a 
Fermi gas has been discussed in Ref . [45| , but in this case 
the backbending is expected to take place at \k\ 3> k F . 
The BCS-like dispersion, suppression of DOS near zero 
energy, and backbending of the dispersion suggest that 
the gas is in the pseudogap regime for 77 = —0.5 and 
77 = at least up to T = T F . 

For rj = 1.0, we find that the T-matrix does not 
have any poles (for T/T F = 0.5 and 1.0) as the bound 
state is pushed against the two-particle continuum in 
the molecule spectral function. From Fig. |H we observe 
that the spectral function (for 77 = 1.0) at large mo- 
menta corresponds to free-particle-like excitations. The 
low-momentum part is more interesting since it is best 
described in terms of non-interacting polarons. In this 



regime (|fe| < k F ), the spectral weight is centered around 
energies that coincide with the dispersion of the attrac- 
tive polaron. We define the (attractive) polaron energy 
as in Refs. [2l|, |4f| |47[ , that is, we consider a single spin-| 
impurity embedded to a Fermi sea of spin-'f atoms. The 
polaron dispersion follows from equations 

w(fc) + ^-e fc -ReS 4 (fc J w(fc))=0, (20) 
u p {k) = w(fc) + /4, (21) 

where fj,\ is tuned such that state | \) is empty. To en- 
sure that the calculation is self-consistent, we set fif = (x, 
where /j is solved from the number equation (|14p corre- 
sponding to the original balanced system. Equations {20 1 
and (|21[) generalize the zero-temperature analysis [21 1 
to finite temperatures. Moreover, at T = 0, Eqs. ([20]) 
and (|2ip are equivalent to the analysis based on a varia- 
tional wave function for the polaron [48|, [49| . 

Our analysis suggests that outside the pseudogap 
regime, even a balanced 2D Fermi gas can be effectively 
described as a gas of nonintcracting polarons. For 2D 
systems, the polaron description has been used by Zhang 
et at ,4j as a possible explanation for their experimen- 
tal observations. For 3D systems, similar speculations 
have appeared earlier in Ref. [50| . Our findings thus sup- 
port the scenario put forward in Ref. |^| and indicate a 
crossover (at fixed temperature) from a polaron gas for 
small attraction to a pseudogap regime at strong attrac- 
tion. 



B. Final state — impurity 

The quasiparticle excitations of the final state have 
been thoroughly discussed in Ref. [21( at T = 0. The 
main features in the final-state spectral function corre- 
spond to the attractive and the repulsive polaron |2l], |47|, 
l5l| . while the contribution from a bound state carries 
only a small spectral weight. Here we briefly discuss the 
qualitative changes that the finite temperature induces. 



As in Rcf. [21], we first assume that the initial mixture 
of spin-f and spin- J, atoms is noninteracting. The chem- 
ical potential is then given by the ideal gas expression 
H (T) = k B T \n(e £ F /k B T - l). The finite temperature 
results in two qualitative changes: the attractive polaron 
acquires a finite lifetime and the molecule-hole continuum 
merges with the attractive polaron branch. Furthermore, 
as the temperature increases, the spectral weight is dis- 
tributed more equally between attractive and repulsive 
polarons. We note that the properties of 2D polarons 
have recently been studied experimentally in Refs. @, Q 
and the experimental data agrees reasonably well with 
the non-self-consistent T-matrix calculation [3, ]2m . 
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FIG. 6. (Color online) Spectral function of the final-state 
atom (impurity) when the majority atoms are dressed by 
strong initial-state interactions. The color scale corresponds 
to log 10 Af(k,uj). The strength of the initial-state interac- 
tion is denoted by »7 is in the figure. The temperature is set 
to T — T F and interactions between the final-state atom and 
spin-t atoms correspond to rj {s = 0. The level splitting be- 
tween states | f) and |/) is denoted by A and is assumed to 
be positive. 



ifested as a suppression of the spectral weight associated 
with the polaron states. 



V. RADIO-FREQUENCY SPECTROSCOPY 

We consider radio-frequency spectroscopy in which a 
balanced mixture of atoms in hypcrhnc states 1^) = | ^) 
and \i 2 ) = | \) is coupled to rf photons inducing a tran- 
sition | \) — > |/) j52[. Within the linear response frame- 
work, one obtains 

4f (Wrf) = -TT^rf Imx fl (fc = 0, -w rf - m + M/)- (22) 

At finite temperatures, the retarded correlation func- 
tion Xr can be computed from the corresponding time- 
ordered correlation function as [531 ] 

Im Xi?( fc - w ) =tanh(-/3w)Imx(fe,o;), (23) 

where 

X(r-r',t-t') = 

-i9{t-t'){[^){r,t)^{r,t)^l{r',t')^ f {r',t')]). (24) 

Since we evaluate the correlation function for w = — w rf — 
IX + fif <C — e F , we have tanh(i/3w) « —1. We neglect 
the vertex corrections I54|-[5rl and obtain 



X (k=0,in n ) 

1 



E 



dg 



Gf(Q^^ m )Gi(q,i^ n + i^ m )- (25) 



Using the spectral representation for the Green's func- 
tions, the analytic continuation can be performed exactly 
and we obtain 



JrfH 



|tt rf l 2 



dg 

(2tt) 2 



dz 
2^ 



A i(q,z 



n/ 



w.) 



x A f (q, z) [n F (z -u- A ff - m) - n F (z)] , 



(26) 



Next we consider the case where the impurity (final- 
state atom) is dynamically created from an interacting, 
balanced gas of spin-f and spin- 4- atoius. The chemi- 
cal potential is now determined by the interacting initial 
state and it is solved from the number equation (fi"4|). In 
Fig. [6l we show the spectral function of the final state 
atoms for different initial-state interactions. For illustra- 
tive purposes, we fix the final-state interaction to r] is = 0. 
In contrast to the experiments in Refs. [J, la], this final- 
state interaction is in the strongly interacting regime. 
With increasing initial-state attraction, the initial state 
becomes increasingly paired, and the added impurity is 
less likely to distort the cloud of majority atoms. Thus 
the impurity behaves more and more like a free particle 
as one approaches the BEC limit. In Fig. [BJ this is man- 



where A^f is the level splitting between states | \) and 
|/) and we have defined the detuning u = w rf — A+f- 
If interactions between states | f) and |/) are negligible, 
then the final-state spectral function becomes AJq, z) = 
2tt5(z — e q — A^ f ). Here we have set /i, = since the final 
state is initially empty. We have also taken explicitly into 
account the level splitting A^. The rf current takes the 
form 



J rf (w) 



TTl^fl 2 



j^ 4i.(g. e „-H-u) Me«-^-«), 

(27) 
where we have neglected the term n F (s + A^*) corre- 
sponding to the occupation of the final state. At finite 
temperatures this term is, in principle, nonzero but we 



which renders np(e^ + A 



l 17' 



ncg- 



assume A^* ^> e F , 

ligiblc. In the current experiments [J, |6j], the Zeeman 
splitting between the relevant hyperfine spin states is 
A +f ~ 50 MHz 57 1 , whereas the Fermi energy is in the 
250 kHz [J, Q . Hence, we expect that 



range e F ~ 10 



neglecting the occupation of the final state is justified 
even at temperatures T ~ T F . 

We note that the more general form of the rf current in 
Eq. ([26]) is phenomcnological since all vertex corrections 
are excluded. Therefore our formalism does not include 
bound-to-bound transitions [25|, [55( ■ Since Eq. (|2"6) con- 
tains information about the final-state polarons, we will 
use it to analyze the NCSU experiment. 



A. Quasi-2D geometry and experimental 
parameters 

To compare our results to the experiments [4j, |6[ , we 
need to establish a connection between the dimension- 
less interaction parameter 77 and an external magnetic 
field B which is used to tune the interactions via the Fes- 
hbach resonance. Since the experiments are performed in 
a quasi-2D geometry, we calculate the two-body binding- 
energy e B in terms of the 3D scatte ring length a 3D using 
the quasi-2D form of the T-matrix |58l - l6C| , 



75 M 



2v27r/m 



4/%D - H- U / U z) ' 

where function T(z) is given by 
du 1 



F(z) = 



4^u 3 / 2 



ko-- 



-2u s 



(28) 



(29) 



The trap frequency along the tightly confined direction 
is denoted by w z and the length scale associated with the 
confinement is given by £ z = \J\/muJ z for atoms with 
mass to. The two-body binding energy e B corresponds 
to a pole of the vacuum T-matrix and is determined by 
the equation (we assume e B > 0) 



£Ja 3D -T(e B /u Jz ) = 0. 



(30) 



We solve Eq. ([30]) numerically to obtain e B = 
£ B(4/ a 3i))- To convert the external magnetic field to 
the 3D scattering length a 3D , we use the parameters of 
Ref. [53] • The two-body binding energy is converted to 
the dimensionless interaction parameter 77 using Eq. (|16[) . 



B. MIT experiment @] 

The MIT experiment (y) is performed using the three 
lowest sublevelsll), |2), and |3) on the 6 Li ground- 
state manifold [57j. The set of initial states is given by 
{I t)> I I)} = {| I)? |3)}i an d the final state corresponds 
to |/) = J 2). Interactions between atoms in the two 



initial states are enhanced by a Feshbach resonance at 
B = 690.4 G. The initial state | \) coupled with the 
final state |/) is chosen such that the final-state inter- 
actions are minimized [6J. In the MIT experiment, the 
two-body binding energy of the final state dimer is much 
larger than the initial-state dimer energy, e' B 3> e B . As 
in the 3D case J25|, the final-state interactions are im- 
portant only if rj^f corresponding to e' B is close to the 
2D unitarity r\ = 0. The MIT experiment is outside this 
regime and hence we neglect the final-state interactions 
altogether. We note that the final-state interactions are 
important for the bound-to-bound transition in the rf 
spectrum [29j . but we do not consider this transition in 
the present work. Although the MIT experiment stud- 
ied the whole dimensional crossover from 3D to 2D, we 
concentrate on the 2D limit of the experiment. Using 
the methods of Ref. [6CJ , it is possible to study rf spec- 
troscopy across the whole dimensional crossover. Due to 
the multiple bands involved with this formalism, the cal- 
culation becomes numerically more demanding and we 
postpone such studies for future research. 



B[G] 



690.7 
720.7 



'/\. 



-0.504 
-0.364 



TABLE II. Interaction parameter r\ corresponding to the 2D 
limit of the MIT experiment. The values of the experimental 
parameters can be found from Ref. JfJ. To convert the exter- 
nal magnetic field to the 3D scattering length a 3D , we have 
used the parameters from Ref. [571 ]. 

The dimensionless interaction parameter 7/ for the 
initial-state interactions is shown in Table iHl for the rele- 
vant values of the external magnetic field. The tempera- 
ture at the MIT experiment was estimated to be of the or- 
der of T F and the peak Fermi energy was approximately 
s F = h x 10 kHz p ■ We note that the data corresponding 
to the deepest optical lattice (the regime of our interest) 
was deeply in the 2D regime since e F /hjj z sa 0.04. Thus 
we expect that this limit is well described by our 2D cal- 
culation. 

Since the final-state interactions mainly affect the tails 
of the rf spectra and round off the onset of the pairing 
peak [0, [29j, the qualitative features of the rf spectra 
measured in Ref. p can be understood by studying the 
spectral function of the initial state. In the absence of 
final-state interactions, Eq. ([27[) suggests that the broad 
peak ("pairing peak") at positive detuning arises from 
the incoherent part of the spectral function correspond- 
ing to the dissociation of the preformed pairs (see Fig. @|. 
In Sec. lIV Al we showed that the threshold for the incoher- 



ent part is approximately given by ui th (q) 
Using Eq. ([27]) we obtain a threshold 



ir(<7)+M- 



-w t h(9 = °) - M 



■ a (q = 0) - 2/j. (31) 



Thus the onset of the pairing peak is related to the bind- 
ing energy of the paired fermions. From Fig. we ob- 





FIG. 7. (Color online) Radio-frequency spectra correspond- 
ing to the MIT experiment [g]. The external magnetic field 
is (a) B = 690.7 G and (b) B = 720.7 G. The dimension- 
less interaction parameter r\ for the initial-state interactions 
is given in Table HT1 The temperature is (a) T/T F = 1.5 and 
(b) TjT F = 2.0. The dots indicate the threshold given by 
Eq. (|3ip and the triangles arise from the BCS-like dispersion 
at zero momentum; see Eqs. (|18[) and (1321 1, 



FIG. 8. (Color online) Radio- frequency spectrum at differ- 
ent temperatures corresponding to the 2D limit of Ref. [g]- 
The external magnetic field is (a) B — 690.7 G and (b) 
B = 720.7 G. The dimensionless interaction parameter r\ for 
the initial-state interactions is shown in Table [II] The tem- 
peratures are given by T/T F — 1.0, 1.5, and 2.0 (from top to 
bottom with respect to the peak on the right-hand side) in 
both panels. 



serve that the onset of the incoherent peak is indeed given 
quite accurately by lu . 

In the pscudogap regime, the spectral function is 
peaked around an effective BCS-like dispersion and 
Eq. (|27p can be used to find the corresponding frequen- 
cies in the rf spectrum, 



w ±(<?) =£q±£BCs(9)-M- 



(32) 



However, evaluating co±(q) at q = does not directly 
give the peak location since the spectral function is in- 
tegrated over momentum; see Fig. [7] Furthermore, we 
note that the quasiparticlc dispersion starts to deviate 
from the BCS form for T/T F = 2.0 in Fig.^b). 

The experimentally measured binding energies (the on- 
set of the pairing peak) were found to agree reasonably 
well with the energies of the two-body bound states in 
quasi-2D systems p,|6l[. Our calculation for the chemical 
potential at the regime of the MIT experiment (T > T F 
and j] < 0) predicts [i to be large and negative. Hence, 
the polarization operator Xo S (Qi w ) m Eq. (0 gives a neg- 
ligible contribution to the many-body T-matrix. In this 
limit, the poles of the T-matrix can be solved exactly and 
we obtain e pair (q) 



-£, 



- 2 



/'• 



c B — i,/u, -r 2 C V The onset of the 
pairing peak becomes w (q<) = e B — |e_, which explains 
why the experiment is in good agreement with the two- 
body calculation. In order to observe many-body effects, 
one should either consider lower temperatures or quench 
the system rapidly to the strongly interacting regime as 
suggested in Ref. [601 ]. 

To probe the temperature dependence of the rf spectra, 
we show the calculated rf current at different tempera- 
tures in Fig. [5) We find the best qualitative agreement 
with Ref. Q for temperatures T/T F ~ 1.5 at the Fesh- 
bach resonance (B = 690.7 G) and T/T F ~ 2.0 on the 
BCS side of the resonance (B = 720.7 G). To reach a 
quantitative agreement with the experimental data, the 
rf signal needs to be averaged over the inhomogeneous 
density. Furthermore, one has to use the temperature 



and the Fermi energy as fittingparameters since they 
were not exactly known in Ref. [6j. 

Since the spectral weight carried by the coherent 
branch of the spectral function increases with increas- 
ing temperature (see Fig. 2]), the corresponding peak in 
the rf spectrum at negative detuning becomes stronger 
with increasing temperature (Fig. [8j. Within our the- 
ory, this peak arises from thermally excited quasiparti- 
cles and corresponds to unpaired atoms. For the fairly 
high temperatures required to obtain a qualitative agree- 
ment with the MIT data, the depletion in the density of 
states is qualitatively the same as in Fig. 03a). The gas 
is therefore at the crossover region between the pseudo- 
gap regime and the regime of a normal gas of bosonic 
molecules 1621. 



C. NCSU experiment 01 

In contrast to the MIT experiment, the NCSU exper- 
iment found resonances in the rf spectrum that do not 
correspond to paired fermions but transitions between 
polaronic states. Our analysis in Sec. IIV Al suggests that 
a balanced gas (the initial state) in the weakly attrac- 
tive regime can be considered as a gas of noninteracting 
polarons. Since the atoms in the final state can be de- 
scribed in terms of polarons (see Sec. IIVB[) . the exper- 
imental picture of transitions between polaronic states 
arises naturally. 

We first analyze the line shape of the rf spectrum in 
a homogeneous system for the parameters of the NCSU 
experiment. The dimensionless parameters 77*, and rj^* 
characterizing the initial- and final-state interactions are 
shown in Table Hill We use the values reported in Ref. [J] 
for the initial and the final-state dimer energies and de- 
termine the interaction parameter rj using Eq. (|16p . Since 
the final-state interactions arc weak, wc expect that our 
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B[G] 


Vn 


%f 


719 
810 


-0.330 
0.967 


1.626 
1.767 



TABLE III. Dimensionless interaction parameters corre- 
sponding to the NCSU experiment [J], initial-state interac- 
tions are denoted by J7+1, and 77^ corresponds to the final- 
state interactions. For other magnetic fields and trap depths 
considered in Ref. [J], the initial- and the final-state interac- 
tions are similar and give qualitatively the same results. 



phenomcnological model in Eq. (|26[) is a good description 
of the dressed final-state atoms. 

The rf spectra corresponding to the parameters of 
Table IIIII are shown in Fig. GO For strong attraction 
[Fig. Eta)], the pairing peak at positive detuning uj is 
again characterized by the threshold frequency lj given 
by Eq. ([31~]) . The sharp peak around zero detuning arises 
from the unpaired atoms as in the MIT experiment. We 
find that the inclusion of the final-state interactions tends 
to lower the peak corresponding to the unpaired atoms. 
A crucial difference between Figs. M,&) and (b) is the ab- 
sence of a pole in the T-matrix for the parameters of 
Fig. HJb). Since most of the data presented in Ref. [4[ 
fall to the same regime of initial- and final-state inter- 
actions as Fig. (b), the measured rf spectra cannot be 
explained in terms of confinement-induced molecules. 
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sion of the final-state atoms is effectively replaced by the 
final-state polaron dispersion, and we obtain an estimate 
for the peak location: lo'q = — w ; 
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where uj p {k = 0) is the energy of the final-state po- 
laron. From Fig. [5] (b), we observe that the estimates ui 
and uj'q are indeed quite accurate. 

Similarly to the MIT experiment [(|, the NCSU ex- 
periment has a dual peak structure in the measured rf 
spectra [J| . Since our calculation for a homogeneous sys- 
tem predicts only a single peak, we argue that the sec- 
ond peak is due to the averaging of the rf signal over 
an inhomogeneous density. As in 3D systems |42l444| . 
the auxiliary peak arises from unpaired atoms in the 
low-density region at the edge of the trap. To show 
this explicitly, we average the rf signal over the inho- 
mogeneous density using the local density approxima- 
tion (LDA) [H|. Within the LDA, the local chemical 
potential is given by /j(r) = /x — ^mui^r 2 . The local 
density n LT)A (r) * s calculated using the number equa- 
tion (fTI)) and the density averaged rf signal is obtained as 
I%P A (u>) = /drn LDA (r)J rf (w,r), where / rf (w,r) is the 
rf signal given by Eq. (P5|) for the local chemical potential 
H(r). 

Rather than attempting to exactly reproduce the ex- 
perimentally measured lincshapes, we use the LDA to 
confirm that the auxiliary peak in the rf spectrum near 
the zero detuning is due to trap average. To this end, we 
fix the initial- and the final-state interaction strengths to 
the values given in Table IIIII and take the Fermi energy 

£ f = UJ ± VXot t0 bc 114 - 5 kHz ; see Ref- 0- Tn c total 
number of atoms is given by N tot = Jdrn LDA (r) and 
the peak value of the chemical potential can be tuned to 
reach a desired number of atoms [43|, |63[ • In Fig. [TU1 we 
show the rf current for different numbers of atoms. We 
implicitly assume that the trap frequency is tuned such 
that e f remains constant. 



FIG. 9. (Color online) Radio-frequency spectrum correspond- 
ing to the NCSU experiment [| for (a) B = 719 G and (b) 
B = 810 G (solid lines). The dimensionless interaction pa- 
rameters for the initial- and the final-state interactions are 
shown in Table Hill The temperature is T/T F — 1.5 in both 
panels. The dashed line denotes the rf spectrum in the ab- 
sence of final state interactions. In panel (a), the dot indicates 
the threshold given by Eq. (|3ip . and in panel (b), the dots cor- 
respond to the polaron threshold; see the main text. The sin- 
gle peak structure in panel (b) persists down to T/T F = 0.5, 
which is the lowest temperature considered in this work. 



Based on the analysis in Sec. IIVAI we argue that the 
origin of the single peak in Fig. EJb) is related to the 
effective description of the gas in terms of noninteract- 
ing polarons. In the absence of final state interactions, 
we expect the location of the peak to be given by uj = 

—ojp(k = 0), where u)p (fe) is the initial-state polaron 
energy given by Eq. (|21[) . Furthermore, when the final- 
state interactions are included, the free-particle dispcr- 
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FIG. 10. (Color online) Trap-averaged rf signal corresponding 
to (a) T/T F = 0.5 and (b) T/T F = 1.0. The dimensionless 
interaction parameters are r^ = 0.967 and 7/1 1 — 1.767. The 
rf current is scaled by the total number of particles which is 
given by (a) N tot = 2234, 1791, 1476, 1214, and (b) N tot = 
2073, 1810, 1494, 1190 (from top to bottom). 

On a qualitative level, our theoretical calculation 
agrees with the picture put forward in Ref. [J]: in the 
regime of the experimental parameters, the rf spectrum 
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is best described in terms of transitions between nonin- 
teracting polaronic states. In general, we find that the 
location of the right-hand side peak in the trap-averaged 
rf signal [Fig. [TU] (a)] is not directly related to the po- 
laron energy corresponding to the peak chemical poten- 
tial. Therefore, locally resolved rf spectroscopy [26| is 
desirable to confirm the polaron picture. Within the ap- 
propriate temperature range, our LDA calculation quali- 
tatively reproduces the double-peak structure observed in 
the experiment; see Fig. llOf a). The auxiliary peak due to 
the trap average becomes more pronounced at high tem- 
peratures and renders the polaron peak indistinguishable 
if the total number of atoms is too small [Fig. [TOT b)]. 



VI. DISCUSSION 

In this work, we analyzed pairing in two-dimensional 
Fermi gases above the supcrnuid transition tempera- 
ture. For a gas composed of an equal number of spin- 
f and spin- J, fcrmions, we found evidence of a crossover 
from a nonintcracting gas of polarons to a pseudogap 
regime characterized by a BCS-like dispersion and re- 
duced single-particle density of states near the Fermi en- 
ergy. The details of this crossover as well as the prop- 
erties of the effective polarons in a balanced Fermi gas 
clearly call for further investigations. 

We also analyzed two recent experiments that per- 
formed radio-frequency spectroscopy for 2D Fermi 
gases (jjlg. To take into account the polaronic properties 
of the final-state atoms, we introduced a phenomcnologi- 
cal model that uses dressed Green's functions for both 
initial- and final-state atoms in the calculation of the 
rf spectrum. Although this approach does not include 
the vertex corrections [5J-[5g, our model qualitatively 
explains the different observations in Refs. [J, Q. In par- 



ticular, our polaron-to-pseudogap crossover can be used 
to explain the apparent dichotomy in Refs. [1, Q- 0m 
calculations suggest that Ref. @ probed the pseudogap 
regime, while Ref. [J] provided data from the polaron side 
of the crossover. 

In Ref. |6J, the measured binding energies were found 
to be close to the energies of the corresponding two-body 
bound states. We demonstrated that this effect arises at 
the high-temperature regime, where many-body contri- 
butions to the experimentally measured thresholds in the 
rf spectrum become negligibly small. In order to clearly 
see the effect of many-body corrections to the measured 
binding energies, the temperature should be lower or the 
system should be quickly quenched to the strongly inter- 
acting regime as suggested in Ref. [60| ■ 

Another interesting and theoretically largely unex- 
plored direction is the dimensional crossover from 2D to 
3D. In Ref. [6(, the dimensional crossover was explicitly 
probed, and in Ref. [J], implications of the quasi-2D na- 
ture of the gas seem to be unavoidable since the Fermi en- 
ergy was larger than the level splitting in the tightly con- 
fined direction. The formalism constructed for molecule 
formation in quasi-2D systems [60( has recently been ex- 
tended to study the polaron problem [64j | and could be 
used to investigate the rf spectroscopy of quasi-2D sys- 
tems. 

Note added in proof: Recently, two preprints dis- 
cussing the pseudogap phase in 2D Fermi gases have ap- 
peared |65ll66|. 
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